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THE  EFFECT  OF  FINITE  “BLOB”  SIZE  ON  THE  CURRENT 
CONVECTIVE  INSTABILITY  IN  THE  AURORAL  IONOSPHERE 

I.  INTRODUCTION 

Ionospheric  irregularities  (or  density  fluctuations)  cause 
scintillation  of  radio  signals  which  has  been  observed  during  diffuse 
auroral  situations  [Fremouw  et  al.,  1977;  Rir.o  et  al.,  1978;  Rir.o  and  Owen, 
1980].  The  irregularities  occur  in  regions  of  soft  particle  precipitation 
(whic.n  are  confined  in  latitude)  and  horizontal  plasma  density  gradients. 
At  present,  the  emerging  picture  is  that  large-scale  density  enhancements 
(known  as  "blobs")  are  produced  in  tne  high  latitude  F  region  and  are 
cor.vected  around  the  polar  ionosphere;  the  sides  of  these  "blobs"  appear  to 
be  the  regions  of  structure  [Tsunoda  and  Vickrey,  1985].  There  has  been  a 
ccnsideraole  theoretical  effort  to  understand  these  irregular ities  in  terms 
of  gradient  driven  fluid  instabilities  [Ossaxow  and  Chaturvedi,  1979:  Huba 
and  Cosa,<ow,  1?8C;  Vickrey  et  al.,  1980;  Keskinen  et  al . ,  '?80:  Chaturvedi 
and  Ossakow,  1981,  1983;  Keskinen  and  Ossakow,  1983;  Satyar.arayar.a  and 
Ossakow,  1983;  Gary,  1 98 U ;  Huba,  1984;  Satyanarayar.a  et  al.,  1985].  Ir. 
particular,  the  current  convective  instability  [Lehr.ert,  1958;  Kadomtsev 
and  Nedcspasov,  I960]  has  been  suggested  as  a  generation  mechanism  of 
plasma  irregularities  in  the  high  latitude  ionosphere  [Ossakow  and 
Chaturvedi,  1979:  Fejer  and  Kelley,  1930;  Har.uise  et  al . ,  1931;  Vickrey  and 
Kelley,  1983].  This  instability  can  be  excited  in  a  weakly  collisior.al, 
magnetized  plasma  which  contains  a  field-aligned  current  and  a  transverse 
density  gradient.  This  ccnf iguratior.  exists,  at  times,  in  the  auroral 
ior.ospnere  along  the  sides  of  density  enhancements  ("blobs"). 
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The  theoretical  study  of  the  current  convective  instability,  as 
applied  to  the  auroral  F  region,  has  evolved  over  the  years.  Initial 
studies  usee  the  local  approximation  which  is  valid  for  modes  such  that  kL 

>>  1  where  k  is  the  wavenumber  and  L  is  the  scale  length  of  the  density 
gradient  [Ossakow  ar.d  Chaturvedi,  1979;  Vickrey  et  al.,  1980;  Chaturvedi 

ar.d  Ossakow,  1931;  Gary,  1984].  Subsequent  studies  have  considered  a  host 
of  nonlocal  effects  which  can  be  important  and  are  neglected  in  a  local 
analysis,  e.g.,  magnetic  shear  [Huba  and  Ossakow,  1980];  velocity  shear 
[Satyar.arayana  and  Cssakcw,  1983];  long  wavelength  modes  (i.e.,  kL  <  1) 
[Huba,  1984];  finite  current  channel  width  [Satyar.arayana  et  al.,  1985]. 
Overall,  these  effects  have  a  tendency  to  reduce  the  growth  rate  of  the 
instability,  with  the  exception  of  velocity  shear  which  can  stabilize  short 
wavelength  modes  (these  with  kL  >>  1).  Despite  the  modest  reductions  in 
growth  rate  the  current  convective  instability  has  remained  a  viable 
mechanism  to  generate  density  irregularities  with  scale  lengths  1-10  km. 

A  potentially  important  nonlocal  effect  that  has  not  been  considered 
to  date  is  the  finite  extent  of  toe  "blob”  along  the  ambient  geomagnetic 
field.  For  example,  recent  studies  of  barium  aloud  dynamics  have  indicated 
that  the  parallel  extent  of  the  cloud  along  the  field  can  affect  its 
stability  properties  [Goldman  et  al.,  1976;  Sperling  et  al.,  1984;  Sperling 
ar.d  Glas3mar.,  1985;  Drake  et  al.,  1985].  Physically,  the  reason  that  this 
effect  could  be  significant  is  the  following.  The  current  convective 
instability  requires  density  and  potential  fluctuations  both  parallel  and 
perpendicular  to  3Qt  i.e.,  assuming  a  Fourier  expansion  of  the  modes  we 

have  <  »  k  e  *  k  e  where  k  is  the  wavenumber.  Based  on  local  theory, 
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the  instability  attair.3  maximum  growth  when  k  /k  -  io  /o  where 
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j  and  c  are  the  plasma  conductivities  parallel  ar.d  perpendicular  to  the 

I  1 


|"3 


O 


field,  respectively.  For  typical  F  region  parameters  we  note  that 

8  4 

a  /a  -  10  so  that  A  -  10  A  where  A  is  the  wavelength  of  the  mode. 

II  II 

Thus,  for  fluctuations  with  A  -  1-10  km,  local  theory  requires  parallel 
wavelengths  such  that  A  ^  -  1QJ  -  IQ3  km  for  maximum  growth.  On  the  other 
hand,  the  parallel  extent  of  a  "blob"  is  typically  only  several  hundred  km 
so  that  local  theory  is  probably  inadequate  to  properly  describe  the 
current  convective  instability  in  the  auroral  ionosphere  (at  least  for  the 
fastest  growing  modes  of  interest). 

In  this  paper  we  develop  a  nonlocal  theory  of  the  current  convective 
instability  which  considers  the  finite  extent  of  an  ionospheric  "blob" 
parallel  to  the  geomagnetic  field.  We  show  that  the  physical  picture 
presented  in  the  previous  paragraph  is,  in  fact,  reasonably  accurate.  For 
a  mode  that  has  a  parallel  structure  sufficiently  short  to  "fit"  the 
fastest  growing  mode  into  the  "blob" ,  we  recover  the  maximum  growth  rate 
predicted  by  local  cheery.  However,  modes  that  have  a  parallel  structure 


determined  bv  the  extent  of 

the 

"blob"  along 

the  field  oar.  have  t 

:rowth 

rates  substantially  -educed 

from 

the  maximum 

growth  -ate  expected 

from 

local  theory.  For  typical  auroral  ionosphere  parameters,  the  reduction  in 
the  maximum  growth  rate  of  the  instability  for  medium  scale  irregular ities 
(1-10  km)  can  be  one  to  two  orders  of  magnitude. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
present  the  physical  model  and  basic  equations  used  in  the  analysis.  In 
Section  III  we  derive  the  dispersion  equation  for  the  current  convective 
instability  which  explicitly  includes  the  finite  extent  of  the  "blob"  along 
the  ambient  magnetic  field.  In  Section  IV  we  present  analytical  and 
numerical  results.  Finally,  in  Section  V  we  summarize  our  findings  and 
apply  our  results  to  the  auroral  ionosphere. 
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II.  MODEL  AND  FUNDAMENTAL  EQUATIONS 

We  consider  a  slab  geometry  and  a  plasma  configuration  as  3hcwr.  in 

Fig.  1  .  The  ambient  magnetic  field  and  current  are  in  the  z-directior. 

(3  -  Bn  e  and  J  -  J  e  ).  The  ambient  density  profile  used  is  given  by 
"*  U  Z  *"  u  z 


nb(x) 

0  <  z  <  ZQ 

n 

“i 

otherwise 

where  nb  denotes  the  "blob"  density  which  is  confined  to  a  limited  region 

along  3Q  and  is  inhomogeneous  transverse  to  3^ ,  and  n<_  denotes  the 

homogeneous  background  ionosphere.  We  consider  low  frequency  fluctuations 

such  that  3/3t  <<  2  ,  v  where  2  and  v  represent  the  gyrof requercy  and 

ci  an  a  ctn 

collision  frequency  with  neutrals  associated  with  the  a  species.  We  also 

assume  v  /a  «  \> .  /a.  «  1  which  is  relevant  to  the  F  region  of  the  ior.o- 

spr.ere.  Finally,  for  simplicity  we  assume  a  cold  plasma  (i.e.,  Ta  -  T^  * 

and  that  the  ambient  current  is  carried  oy  electrons  (i.e.,  =  -  na 

Vn  e  ). 

Oz  z 

Within  the  context  of  our  assumptions,  the  basic  equations  of  our 


analysis  are  continuity,  momentum  transfer,  charge  neutrality,  and  Ampere's 


(5) 


V-J  -  7.[r.e(yi  -  yj]  -  0 


7x3 


4tt 
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J 


(5) 


where  the  variables  have  their  usual  meanings  and  we  are  ir.  the  neutral 

frame  of  reference  (i.e.,  V  -  0).  We  take  the  electric  and  magnetic 

-n 

fields  to  be  represented  by  potentials  as 
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1  3Az  ^ 

-  7*  -  c  ST  e; 


(7) 


and 


B 


7A  x  e 


(3) 


where  B0  is  the  ambient  field,  and  $  and  A2  are  the  electrostatic  and 

vector  oocentials,  -esoecti vely .  We  consider  only  A,  since  J  >>  J  . 
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The  electron  cross-field  motion  is  given  by 


v 

-n 


!y 


x  e 


while  the  parallel  motion  is  given  by 


v 

91 


e  en 


,  3  a 

tVssri- 


(9) 


(10) 


The  ion  cross-field  motion  is  given  by 


A  V  . 

G  i  n  <-» 

v,  »  -  r-  7  <|>  x  e  -  ^  V  i 

-i±  3  _l  2  ai  3  1 


and  che  parallel  motion  is  giver,  by 


a  v 

v.  -  -  —  — —  V  <<  V 

H  ai  vin  e !  el 


We  row  substitute  (9)-('2)  into  (2),  (5),  and  (5)  and  arrive  at  the 
fallowing  equations: 


-  |  7  f  x  e  *7n  ♦  r —  |-  72A  -  0 

3t  B  j/  z  4ire  3z  ^  2 


v . 

i  n  ^ 


-=£  |  7  •  ( r.7  4)]  +  72A  -  0 

Q.  3  v  r1  4-ire  3z  i  z 


4ir  +  I 
cn  3z  c  3t 


’A^a  f|  m  ,71  ;  /  r  t?  . 

3  " 'hi  an 


The  equ:.l  tori  urn  solution  of  (13)~0  5)  requires  that 


c  ~  n  5  2 

-  £  7  h  x  e  •  7-  «•  r-rr  h  ?  A_„  »  0 


3  j_  0  z  "0  4ire  3z  j_  z0 


v . 

L  G  G  r 

~  f  7  •  r.,7  ® 

u .  3  0  i  J 


—  I-  72A 
Ut9  jz  j_  zC 


v2a  -  J*i-  I!° 

j_  zO  cna  3z 


We  have  assumed  r.3  -  n^'x.z)  [see  U)j  and  will  also  ta*e  and  Az0  to  be 
functions  or._y  of  x  and  z.  [Neglect  of  the  y  dependence  or.  $  is 


equivalent  to  assuming  -n  -  0;  we  note  that  it  is  the  y  component 
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of  Eq  which  car.  drive  the  usual  E  x  B  gradient  drift  instability  in  the  F 
region  ar.d  we  are  excluding  this  possibility.]  From  (16)  and  07)  we  then 
find  that 


c  l_ 
Uire  3z 


V^A 
Vj_  zO 


0 


(19) 


This  potential  generates  an  equilibrium  electric  field 


“3xni/nb( 


X;  e. 


0  <  z  <  z0 


E'  e  otherwise 

Oz  z 


(2«) 


which  causes 

inhomogeneous 

"blob”  plasma 

drift  V'  -  - 
-Oz 


a  parallel  electron  drift  V  =»  -  ( eE^/m^v^ )e 


0  z  e  er.  z 


cross-field  drift  V  (x)  »  -  [cE.  /B)(n,/r.  . (x)]e 


~0j_ 


OX 


In  the  background  plasma  there  is  only  a  parallel 

f eE '  /mu  ]e  . 

'  Oz  e  er.  z 


and  an 
in  the 
electron 


III.  DISPERSION  EQUATION 

We  linearize  C 1 3 ;  —  { : 5 )  in  order  to  obtain  the  dispersion  equation 
which  describes  the  current  convective  instability.  We  use  the  equilibrium 
derived  Lr.  Section  II  ar.:  assume  perturbed  quantities  vary  as  p  *  p(x)  exp 
' Yt  -  The  linear ; zed  equations  are 


o  -  A  £  *  -  £  3  * 

-  £  7  $  x  a  *7  r.  -  -  7  0 _  x  e  *  Vn  -  — -  m-  7  A 

3t  31  zj_0  3  ^  0  z  5:  j_  z 


(25) 


~  3  7‘ 


'■oV 


£  7 
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0' 


_c  __3  7  2- 
4ire  3z  j_'  t 


'  26 ' 


and 


—  3  <1 

72  A  -  iii®  :v  *  —  —  -  -  — 
■z  c  Oz  cn  3z  c  3t 
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(27) 


In  tne  following  we  will 
for  <  7.  <<  1 V,  3/3z 1 . 


ignore  terms  proportional  to 
Retaining  these  terms  leads  t 


V  i.  which  is 
1  3 

0  two  effects: 


valid 

(1  ) 


a  real  frequency  is  produced  proportional  to  where  x  -  xQ  is  the 
location  of  mode  localization,  and  (2)  velocity  shear  stabilization  cf 
short  wavelength  modes  (Satyanarayar.a  and  Ossakow,  1983).  We  are 
intentionally  ignoring  these  effects  in  order  to  highlight  the  influence  of 
finite  "blob"  size  on  the  current  convective  instability. 

Eliminating  n  from  these  equations  we  obtain  two  coupled  differential 
equations  for  A^  and  5, 


[y 


(28) 


and 


n  i!i* 
a  Dr  c  3z 


(29) 


where  D„  -  (c  /w  jv  is  the  resistive  diffusion  coefficient, 
pe  e 

uPe  “  ^noe2/me’  a  M  (Qe/ve^ai/vi rj»  and  Y0  “  “(ve/fle^  V0z  3Zr 1  no/3x* 

In  writing  (28)  we  have  taken  (v  /y)0/3z  <<  1  which  can  be  shown  a 

u  z 

posteriori;  this  term  dees  not  affect  the  growth  rate  of  the  mode  out  only 
causes  the  mode  to  have  a  real  frequency.  Finally,  combining  (28)  and 
(29),  we  obtain  the  mode  equation  for  the  current  convective  instability  in 
terms  of  the  vector  potential  A^, 


3_ 

3z 


3A  Y  _  3A 

Dr  -  1  r  *  °r\  if  '  <Y 


(30) 


Prior  to  solving  (30)  for  the  density  profile  shown  in  Fig.  lb,  we 

first  consider  a  blob  of  infinite  extent  along  the  field.  We  Fourier 

expand  modes  parallel  to  B0,  i.e.,  p(z)  -  p  exp(ik  z)  and  use  the  local 

^  z 

approximation  [k  (3  In  n  / 3x)  1  >>  l]  so 


that  (30)  can  be  solved 


algebraically.  This  allows  comparison  with  previous  results  [Ossakow  and 
Chaturvedi,  1979;  Chaturvedi  and  Ossakow,  1981 j.  The  local  dispersion 
equation  is  given  by 


y 


y 

k  k  /a  D  +  (k* 
Y  y  z  r  } 


0. 


(31 ) 


The  first  term  in  (31)  is  related  to  electromagnetic  effects,  the  second 

term  to  convection  along  the  gradient  (which  causes  the  instability),  and 

the  final  term  to  perpendicular  ion  motion  and  parallel  electron  motion. 

In  the  electrostatic  limit  (k  u  ,  a*  D  >>  y)  the  growth  rate  is  given  by 

y  r  z  r 


Y 


+  Y 


^  (Vky) 

0  ,  2,,  2 
1  +  ax  / k 
2  y 


(32) 


which  agrees  with  the  results  of  Ossakow  and  Chaturvedi  (1979).  This 

growth  rate  has  a  maximum  value  when  k"/k“  =  1/a.  It  is  given  by 

2  y 


Y 

max 


V2 


(33) 


Retention  of  electromagnetic  effects  yields  the  following  growth  rate 


Y  »  -  ^  k2D  +  4  [k2D  *  i»Y_k  k  /a  D  1 
2  r»  -  2  r  Oyz  rJ 


1/2 


(3*0 


*2  2  2 

where  k  -  k  +  ak  and  which  can  be  shown  to  agree  with  Chaturvedi  and 

y  ^ 


Ossakow  (1981). 

We  now  solve  (30)  for  a  density  profile  relevant  to  a  plasma  "blob"  in 


the  auroral  ionosphere.  The  equation  is  solved  subject  to  the  boundary 


condition  -*  0  as  \z\  *  ».  For  a  "top  hat"  density  profile  (Fig.  lb), 
the  solution  to  (30)  in  the  region  z  <  0  and  z  >  zQ  can  be  written  as  a 
plane  wave, 


where 


Az  -  A1  exp(-  1 2 1 ) 


k.  *  k  (— — ) 

1  y  jzz < 


2 

Y  +  k  D 

y  r 

Mote  that  the  subscript  i  denotes  quantities  that  are  evaluated  in  the 
background  ionosphere  'i.e.,  z  <  3  and  z  >  Zq).  The  parameter  R  is  a 
measure  of  the  electrcstatic/electrcmagr.etic  nature  of  the  mode.  For 
Y  <<  k^D  we  note  that  R  -  1  and  the  mode  is  essentially  electrostatic.  In 

y  r 

the  opposite  limit,  Y  >>  k^O^  the  mode  is  essentially  electromagnetic  with 
R  <<  1  . 

In  the  region  0  <  z  <  zQ,  i.e.,  region  occupied  by  the  density 
■enhancement ,  the  solution  to  A„  is  given  by 


A2  -  A^  exp( ikj z)  *  A,*  exp(ik2z) 


where 


■■>v  »>' 


V  *.  A  *.  \ 


Again,  the  subscript  b  denotes  quantities  that  are  to  be  evaluated  in  the 


blob  (0  <  z  <  Zq). 

To  determine  the  dispersion  equation  we  match  the  plane  wave  solutions 

at  z  =  0  and  z  -  Zq.  The  appropriate  matching  conditions  are  obtained  from 

(23)  and  (29).  Namely,  we  require  that  A  and  aD  3A  /3z  be  ccntir.ous 

z  r  z 

across  each  boundary.  Using  these  boundary  conditions,  we  find  that  the 
dispersion  equation  is  given  by 


where 


exp(iikz3)  *  A 


1  Mr  r= 

_ 

k1 

i  ♦  ir  7- 
k. 
i 


1 

ir  7~ 
k . 
l 


1  - 


*2 

ir~ 

l 


Ak  -  k 


-L  rls  -}uz 


) 


(42) 


(u3) 


(44) 


and  kj_,  k1  ,  and  k2  are  given  in  (36),  (39),  and  (40),  respectively. 


IV.  ANALYTICAL  AND  NUMERICAL  RESULTS 

We  now  present  analytical  and  numerical  results  for  the  growth  rate  of 
the  current  convective  instability  based  upcn  (41).  We  cast  (41 )  into 
dimensionless  form  and  obtain 


exp( ikyZgAk) 


(45) 


where  k^  *  ky^r>  zq  *  Zg^/a  “  Ak /k  and  Lr  -  Dp/Y^  is  the  resistive 

diffusion  length.  First,  we  note  that  in  the  short  wavelength  limit,  i.e., 

A  A  A 

when  kyZg  >>  1 ,  which  implies  Ak  -  0,  so  that  A  -  1  i3  a  solution 

to  (45).  The  growth  rate  is  given  by 


Y 


(46) 


In  the  electrostatic  limit  (y  <<  k2D  }  we  note  that  R  -  1  and  Y  -  Y ./2 

y  ry  0 


which  corresponds  to  one  maximum  growth  rate  of  the  instability  based  upon 
local  theory  I  see  (33)2.  This  result  makes  sense  physically  since  the 


limit  kyZg  >>  1  corresponds  to  modes  with  parallel  wavelengths  much  smaller 


than  the  parallel  length  of  the  blob;  hence,  the  outer  regions  (i.e.,  z  <  0 
and  z  >  Zg)  have  little  affect  on  the  dynamics  of  the  instability. 


On  the  other  hand,  for  long  wavelengths,  i.e.,  k  zn  <<  1,  tne  growth 

y  u 


rate  is  3trongly  affected  by  the  finite  length  of  the  blob.  The  growth 
rate  is  reduced  and  k2  >>  k1  .  In  this  limit  A  -  -  1  and  the  dispersion 
equation  becomes 


k  z„  Ak  -  (2m  +  1  )ir 


(47) 


where  a  *  I,  2,  ...  is  an  integer.  Equation  (47)  yields  the  following 


We  note  that  this  growth  rate  is  independent  of  the  value  of  R.  Also,  the 
fastest  growing  mode  in  this  regime  corresponds  to  m  -  0. 

In  Fig.  2  we  present  numerical  results  based  upon  (45).  'We  plot 

A  A  A 

Y  -  Y/Yq  versus  k,  -  k.^L^  for  zQ  -  Zg//a  -  1  and  r  ■  0.2  for  several 

values  of  mode  number  m.  The  solid  curves  are  numerical  results  while  the 
dashed  curves  are  based  upon  the  analytical  results  (46)  and  (43).  Note 
the  excellent  agreement  between  these  results  in  the  short  wavelength 

A  A 

(i.e.,  ky  >>  2*)  and  long  wavelength  (i.e.,  k^  <<  2ir)  limits.  In  the 
short  wavelength  regime  the  modes  are  electrostatic  (R  -  1)  and  have  a 
growth  rate  which  is  the  same  as  the  maximum  value  obtained  from  local 


theory  [y  -  7-/2,  see  (33)].  In  the  long  wavele.ogtn  limit  the  modes  are 

vj 

electromagnetic  (?.  <<  1)  and  have  a  smaller  growth  rate  than  in  the  sr.ort 


wavelength  limit. 


Physically,  the  reduction  in  growth  car.  be  explained  as  follows.  The 

current  convective  instability  requires  a  density  and  potential 

perturbation  parallel  to  3Q  which  corresponds  to  a  parallel  wavenumber  in 

the  local  analysis  (the  second  term  in  (31)  is  the  driver  and  depends  upon 

Yq  and  k2).  From  local  theory,  the  instability  achieves  maximum  growth 

when  k_/k„  »  l//a  *  fa  /a  }  where  a  and  a  are  the  oerper.dioular  and 
z  y  'll  1  I 

«  Q 

parallel  conductivities,  resoectively .  In  the  ior.osohere,  o  /a  ■  10  « 

’  'll 

1  so  that  the  instability  favors  very  long  wavelengths  parallel 


to  3q  relative  to  the  perpendicular  wavelengths  >>  A  ) .  We  define  the 
maximum  effective  parallel  wavenumber  that  car.  be  supported  within  the  blob 


kz(eff) 


*  2ir/z 


0* 


(49) 


Using  (49)  to  determine  the  value  of  k^zQ  for  maximum  growth  based  on  local 
theory  /k^  »  /a]  we  find  that  k^zQ  »  2tt.  Thus,  for  k^zQ  >>  2m  we  expect 
that  the  blob  is  sufficiently  lcr.g  to  support  a  parallel  wavelength  which 
yields  maximum  growth,  i.e.,  k,  >  !<2(eff)-  This  is,  in  fact,  quite 

A  A  A 

apparent  in  Fig.  2  where  Y  -  0.5  for  k^z^  >>  2ir.  However,  for  modes  such 

A  A 

that  ky2g  <<  2ir  it  is  not  possible  for  any  mode  to  satisfy  the  condition 

k  /k  -  T/Vot.  Therefore,  the  modes  grow  at  a  reduced  growth  rate;  acain 
*•  / 

this  is  apparent  from  Fig.  2.  In  fact,  this  corresponds  to  the  limit 

/a  k  /k  >>  1  based  on  local  theory.  Making  use  of  this  approximation  in 
z  y 

(32)  we  find  that 


Y 


(50) 


which  agrees  with  (48)  if  we  make  the  identification  k 


-1 


,/  (2m  +  )ir. 


V.  DISCUSSION 

We  have  studied  the  effects  of  finite  size  of  high  latitude  "blobs" 
parallel  to  magnetic  field  on  the  current  convective  instability  in  auroral 
ionosphere.  We  find  that  the  growth  rate  of  the  instability  can  be 
substantially  reduced  in  the  case  of  a  finite  size  "blob"  from  the  value 


obtained  when  the  "blobs" 


are  assumed  to  be  infinite. 


For  short 


wavelengths,  such  that  kyZg  >>  /a  (where  k  is  the  wavenumber,  Zg  the  blob 

size  ar.d  a  -  Q  Q./\>  v.  -  o  /a  ),  the  growth  rate  remains  relatively 
M  S  1  '  1  II 

unaffected  (especially  in  the  electrostatic  case)  by  the  bloc  size.  In  the 
long  wavelength  case,  when  kyZg  <<  /a,  the  growth  rate  becomes  proportional 
to  the  wavenumber  and  tne  blob  size.  Thus,  for  a  given  blob  size,  the 
longer  transverse  wavelengths  have  smaller  growth  rates  and,  similarly,  for 
shorter  blob  sizes,  the  growth  rates  are  reduced  [3ee  (13) j. 

The  current  convective  instability  has  been  discussed  in  the  auroral 
ionosphere  recently,  where  field-aligr.ed  currents  are  a  constant  feature 
and  plasma  density  enhancements  ("blobs")  have  been  observed  with 
structured  walls.  We  now  apply  our  results  for  the  instability  to  this 


situation.  The  typical  parameters  in  this  situation  are  the  following: 
ambient  density  r.^  -  *0®  cm-^,  bloc  density  •%  -ox  *0°  cm-3,  field- 
aligned  current  density  J  ,,  <  102  uamp/m2  (which  corresponds  to  an 
electron  oarallel  drift  v„  <  500  m/sec  for  r.«  -  10°  cm-21),  the  scale 

W  |  -  W 

length  associated  with  tne  "bloc"  gradient  lrI  -  50  kms,  the  ion-neutral 

,4 


:ili3icn  frequency  v . ^  -  5  x 


-  O  -  * 

~  3<2C 


v.^ccron  c  c  _ :  3  i  cr.  frequency 


vo  -  5  x  102  sec  -1  [Tsur.oda  ar.d  Vickrey,  1985].  For  BQ  -  0.5  G,  or.e 

has  v  /Q  -  5  x  10-5  and  v.  /a.  -  1.7  x  10-4  so  that  a  -  fa  a./v  v.  )  - 
e  e  mi  v  e  i  e  m 

10®.  The  maximum  local  growth  rate  is  Y  -  Yft/2  -  3.5  x  10-^  sec-1,  where 

m  0 

Y_  ~  ( v*  /L.J(a  v./v  a.)  -  1.7  x  i0-2  sec-1.  The  "blob"  dimension  parallel 
0"0|N'^eiei' 

to  the  field  is  Zg  -  300  km  CTsur.oda  and  Vickrey,  1985]  ar.d  for  the  above 

set  of  parameters  we  have  <u  -  5.5  x  10'  sec  ,  D_.  -  1.4  x  10  era  /sec, 

pe  "i 

L.  -  D„«/Y,.  -  3  x  i09  cm,  and  z.  -  z  //a  L  -  3.3  x  10-^.  For  A  -  1 
i  ■  i  j  j  u  r  j_ 

km,  k  •  k  L  -  5  x  10®,  ar.d  k  z  -  0.2.  We  find  from  (48),  that  the 
y  y  r  y  o 

r.cr. local  growth  rate  in  this  instance  is  Y  -  0.06  YQ  *  1  .0  x  10-^  sec-1  for 


the  fastest  growing  mode  (m 


0).  This  is  an  order  of  magnitude  smaller 


than  the  maximum  local  growth  Y  ,  i.e.,  Y_  *  8.5  x  10  ^  sec  .  For  a 

m  m 

similar  set  of  parameters,  we  find  that  for  -  10  km,  the  reduction  in 
the  growth  rate  is  ever,  more  substantial,  i.e.,  Y  *  0.006  Y  *1.0  x  10  4 
sec-’.  Therefore,  the  current  convective  instability  growth  rates  may  be 
too  small  for  the  instability  to  explain  the  blob-associated  structure  for 
irregularity  scale-sizes  of  1  -  10  km. 

It  may  be  pointed  cut  that  the  finite  blob-3ize  induced  reduction  in 
the  current  convective  instability  growth  rate  of  the  instability  is  also 
related  to  our  assumption  that  the  instability  region  is  confined  to  the 
"blob"  (i.e.,  the  driving  density  gradient  is  negligible  outside  the  blob, 
namely,  In  the  ambient  ionosphere).  This  is  represented  by  the  requirement 
that  the  modes  decay  exponentially  outside  the  blob.  This  is  a  reasonable 
assumption  in  the  blob-associated  structure  studies.  However,  should  there 
be  a  transverse  gradient  in  the  ambient  ionosphere  at  high  altitude  F- 
region,  and  if  a  positive  correlation  is  observed  between  the  structure  and 
the  field- aligned  currents,  then  the  current  convective  instability  could 
oe  considered  a  cossitie  mechanism  for  scale  sizes  on  the  order  of  a  few  km 
or  less. 
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(b) 

g.  1)  Plasma  configuration  and  slab  geometry.  (a)  Coordinate  system  ar.d 
ambient  plasma  properties.  (b)  "Top  hat"  density  model  for  plasma 
along  the  ambient  magnetic  field. 
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Plot  of  growth  rate  (y  »  Y/Y  }  versus  perpendicular  wavenumber 
(K  *  k^L^}  f>or  3everai  1710(16  numbers  (m  -  0,1,2).  The  solid 
curves  are  numerical  solutions  of  the  dispersion  equation  (45) 
while  the  dasned  curves  are  analytic  solutions. 


use  (46)  while  for  k  <<  l  we  use  (48). 


For  k  >>  1  we 
7 
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0 1 C Y  ATTN  DRSEL-PL-ENV  H.  3CMKE 
0 1 CY  ATTN  J.E.  QUIGLEY 

COMMANDER 

U.S.  ARMY  COMM- ELEC  ENG RG  INSTAL  AGY 
FT.  HUACHUCA,  AZ  35613 

0 1 C Y  ATTN  CCC-EMEO  GEORGE  LANE 

COMMANDER 

U.S.  ARMY  FOREIGN  SCIENCE  1  TECH  CTS 
220  7TH  STREET,  NE 
CHARLOTTESVILLE,  VA  2C9C 
0 ' C Y  ATTN  DRXST-3D 

COMMANDER 

•J.S.  ARMY  MATERIAL  DEV  &  READINESS  CMD 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

O’CY  ATTN  DRCLDC  J.A.  SENDER 

COMMANDER 

•J.S.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  9ACKLICK  ROAD 
3LDG  2073 

SPRINGFIELD,  VA  22150 
0 1 CY  ATTN  LIBRARY 


DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH 
LABORATORY 

A3ERDESN  PROVING  GROUND,  MD  21005 
01CY  ATTN  TECH  LI3RARY, 

EDWARD  3AICY 

COMMANDER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONMOUTH,  NJ  07703 

01 CY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  AL  35309 
01 CY  ATTN  JIM  GAMBLE 

DIRECTOR 

U.S.  ARMY  TRADOC  SYSTEMS  ANALYSIS 
ACTIVITY 

WHITE  SANDS  MISSILE  RANGE,  NM  38002 
0 1 C  Y  ATTN  ATAA-SA 
0 1 C Y  ATTN  TCC/F.  PAYAN  JR. 

0 1 C  Y  ATTN  ATTA-TAC  LTC  J.  HESSE 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
WASHINGTON,  DC  20360 

0 1 C Y  ATTN  NAVALEX  03“  T.  HUGHES 
. 01 CY  ATTN  PME  1 1 7 
0 1 C Y  ATTN  PME  1’7-T 
01 CY  ATTN  CODE  501 1 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CTR 
4301  3'J  ITL  AND  ROAD,  3LDG.  5 

WASHINGTON.  DC  20350 

I*:!  ATTN  MR.  0  Li  3  3  I  N  STIC  12 
01 CY  ATTN  NISC-53 
0 1 C Y  ATTN  CODE  5“04  J.  GALET 

COMMANDER 

NAVAL  OCCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  DA  9  2 1 5  2 

01  CY  ATTN  J.  FERGUSON 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON .  DC  20375 

31CY  ATTN  CODE  >1700  S  .  -  . 
26  CYS  IF  "J NC LASS 
(01CY  I?  CLASS) 

ATTN  CODE  ^730  J.D. 
CYS  I?  UMCLASS,  O'.  C 
01CY  ATTN  CODE  4701  I.  V 


Ossakow, 


.  HUSA,  50 
CY  Ir  CLASS 
VI TKC  VI T3KY 


COMMANDER 

AEROSPACE  DEFENSE  COMMAND/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AF 3 ,  CO  50912 

0  ’  C  Y  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMAND/ XPD 


O'  CY 

ATTN 

CODE 

7500 

DEPARTMENT  OF  THE 

AIR  FORCE 

01  CY 

ATTN 

CODE 

7550 

ENT  AF  3 , 

,  CO  80912 

0  1  CY 

ATTN 

CODE 

7530 

01  CY 

ATTN  XPDQQ 

01  CY 

ATTN 

CODE 

7551 

01  CY 

ATTN  XP 

0  1  C  Y 

ATTN 

CODE 

7555 

01  CY 

ATTN 

CODE 

4730  E.  MCLEAN 

AIR  FORCE  GEOPHYSICS  LABORATORY 

0  1  C  Y 

ATTN 

CODE 

4108 

HANSCOM 

AFB,  MA  01731 

01  CY 

ATTN 

CODE 

4730  3.  RIPIN 

01  CY 

ATTN  OPR 

HAROLD  GARDNER 

20CY 

ATTN 

CODE 

2628 

01  CY 

ATTN  LKB 
KENNETH  S. 

W.  CHAMPION 

COMMANDER 

01  CY 

ATTN  OPR 

ALVA  T.  STAIR 

NAVAL  SPACE  , 

SURVEILLANCE  SYSTEM 

01  CY 

ATTN  PHD 

JURGEN  3UCHAU 

DAHLGREN 

.  VA 

22443 

01  CY 

ATTN  PHD 

JOHN  P.  MULLEN 

01 CY  ATTN  CAPT  J.H.  BURTON 

OFFICER-IN-CHARGE 
NAVAL  SURFACE  WEAPONS  CENTER 
WHITE  OAK,  SILVER  SPRING,  MD  20910 
0 1 C Y  ATTN  CODE  Fjl 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  DC  20375 
01 CY  ATTN  NSP-21 4; 

0 1 C Y  ATTN  NSSP-2722  FRED  WIMBERLY 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
I  AHLGREN  LA  3C  RATS  ?.  Y 
C AHLGRSN ,  VA  22443 

31  C Y  ATTN  CODE  DF-U  R.  BUTLER 

OFFICER  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  222’ 7 
0 1 C Y  ATTN  CODE  435 
C ' C Y  ATTN  CODE  46) 

0 1 C Y  ATTN  CODE  402 
0 1 C Y  ATTN  CODE  420 
01CY  ATTN  CODE  421 


AF  WEAPONS  LABORATORY 
KIRTLAND  AFT,  NM  37117 
01CY  ATTN  SUL 

0 1 C  Y  ATTN  CA  ARTHUR  H.  GUENTHER 
01  CY  ATTN  NTYCE  1 LT .  G.  KRAJEI 

A :  .  A  v 

PATRICK  AF B ,  FL  32925 
0 1 C Y  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AF3 ,  OH  45433 
01 CY  ATTN  AAD  WADE  HUNT 
01 CY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  A  ACC 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  DC  20330 
0 1 C Y  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB ,  MA  01731*5000 
01 CY  ATTN  J.  DEAS 
ESD/SCD-4 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT- PATTERSON  AF 3 ,  OH  45*35 
01 CY  ATTN  NIG  D  LIBRARY 
01CY  ATTN  STD?  3.  3 ALLARD 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AF3  ,  NY  1  344: 

01CY  ATTN  DOC  LI  3 R AR Y / T3L D 
0 ’ 0 Y  ATTN  OCSE  V.  COYNE 

STRATEGIC  AIR  COMMAND/XPFS 
OFFUTT  AF3 ,  N3  63113 

01CY  ATTN  ADWATE  MAJ  3RUCE  SAUER 
0 1 C Y  ATTN  NRT 

01 CY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/SK 
P.O.  BOX  92960 
WCRLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 

0 1 C  Y  ATTN  SKA  (SPACE  COMM  SYSTEMS) 
M.  CLAVIN 

SAMSO/MN 

NORTON  AF3 ,  CA  92409 
( M I NUTEM AN  ) 

01 CY  ATTN  MNNL 

COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HAN3C0M  AFB ,  MA  01131 

01 CY  ATTN  EEP  A.  LORENTGSN 

DEPARTMENT  Of  ENERGY 
LI3RARY  3COM  0-042 
WASHINGTON.  OC  21545 

01  CY  ATTN  00 C  CON  FOR  A.  LA3CWITC 

DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  30X  5400 
AL3UQUSRQUS,  NM  3 " *  15 

01 CY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G ,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  SOX  309 

LOS  ALAMOS,  NM  35544 

0 1 C  Y  ATTN  DOC  CON  FOR  J.  3REEDL0VE 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  303 


VERMORE,  CA 
01 CY  ATTN 

94550 
DOC  CON 

FOR 

TECH 

INFO 

01  CY 

DEPT 

ATTN 

DOC 

CON 

FOR 

L-339 

01  CY 

ATTN 

DOC 

CON 

FCR 

L-31 

R .  H AGE 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  30 X  1663 

LOS  ALAMOS,  NM  87545 


0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

01  CY 

ATTN 

DOC 

CON 

FOR 

R  . 

F.  TASCHEK 

01  CY 

ATTN 

DOC 

CON 

FOR 

E  . 

JONES 

01  CY 

ATTN 

DOC 

CON 

FOR 

J . 

MALIK 

01  CY 

ATTN 

DOC 

CON 

FCR 

R  . 

JEFFRIES 

01  CY 

ATTN 

DOC 

CON 

FOR 

J. 

CINN 

0  1  C  Y 

ATTN 

DOC 

CON 

FCR 

D. 

WESTERVEL 

01  CY 

ATTN 

D  . 

SAPPENF  I 

ELD 

LOS  ALAMOS  NATIONAL  LABORATORY 
MS  D438 

LOS  ALAMOS,  NM  37545 
0 1 C Y  ATTN  S.P.  GARY 
0 1 C Y  ATTN  J .  BOROVSKY 

SANDIA  LABORATORIES 
P.O.  SOX  5300 


ALBUQUERQUE  , 

NM 

571  1  5 

01  CY 

ATTN 

DOC 

CON 

FOR 

W. 

BROWN 

01  CY 

ATTN 

DOC 

CON 

FOR 

A. 

THORNBROUGH 

01  CY 

ATTN 

DOC 

CON 

FOR 

4  • 

WRIGHT 

0’  CY 

ATTN 

DOC 

CON 

FOR 

0  . 

3AHLGREN 

0 1 :  y 

.ATTN 

DOC 

COM 

FOP 

7 1 

4  ' 

0  1  C  Y 

ATTN 

DOC 

CON 

FOR 

SPACE  PRO JSC 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  ?4=50 

0 1 C  Y  ATTN  DOC  COM  FCR  3.  MURPHEY 
CiCY  ATTN  DOC  CON'  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  DC  20545 

01CY  ATTN  DOC  CON  DR.  YO  SONG 


OTHER  GOVERNMENT 


INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO 
AOMIN 

BOL'LOER,  CC  303  0  3 

01  CY  ATTN  0.  C  ROMS  I E 
Q’CY  ATTN  L.  3SRRY 


NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LA3CRATCRISS 
DEPARTMENT  OF  COMMERCE 
30ULDER,  CO  80302 
0 1  C  Y  ATTN  R.  GRUBB 
0  1  C  Y  ATTN  AERONOMY  L A3  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 
P.O.  90  X  92957 
LOS  ANGELES,  CA 
O’ CY  ATTN  I. 

ATTN  T. 

ATTN  V. 

ATTN  S. 

ATTN  D. 


01  GY 
0  1  GY 
01  CY 
01  CY 


90009 

CAR  FUNK EL 

SALMI 

JOSEPHSON 

30WER 

OLSEN 


ANALYTICAL  SYSTEMS  ENGINEERING  CGRP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01SC3 

01 CY  ATTN  RADIO  SCIENCES 


AUSTIN  RESEARCH  ASSOC.,  INC. 
1901  RUTLAND  DRIVE 
AUSTIN,  TX  ?37;8 

0 i C  Y  ATTN  L.  SLOAN 
D '  C  Y  ATTN  R.  THOMPSON 


3ERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  30X  983 
BERKELEY,  CA  94701 


CHARLES  STARK  DRAPER  LABORATORY,  INC 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE ,  MA  02' 39 
01  CY  ATTN'  D.3.  CQX 
01 CY  ATTN  J.P.  GILMORE 


COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  2C734 
01 CY  ATTN  G.  HYDE 


CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

0 1 CY  ATTN  D.T.  FARLEY,  JR. 


ELECTROSPACE  SYSTEMS,  INC. 

BOX  1339 

RICHARDSON,  TX  75080 
0 1 C Y  ATTN  H.  LOGSTON 
0 1 C Y  ATTN  SECURITY  (PAUL  PHILLIPS 


EOS  TECHNOLOGIES,  INC. 

506  Wilshire  Bivd. 

Santa  Monica,  CA  90401 
0 1 C  Y  ATTN  C.3.  GABBARD 
0 1 C Y  ATTN  R.  L2LEVIER 


ESL  ,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 

01  CY  ATTN  J.  R03ERTS 
01 CY  ATTN  JAMES  MARSHAL: 


GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORCE  SPACE  CENTER 
IOC  CARD  3  L/C  KING  CF  PRUSSIA 
P.O.  5CX  5555 
PHILADELPHIA,  PA  19101 
C1CY  ATTN  M.H.  30RTNER 
SPACE  SCI  LAB 


0  1  CY 

ATTN  J. 

.  WORKMAN 

C  *  C  Y 

ATTN  C. 

3  R  ETT I  E 

GENERAL  ELECTRIC  TECH  SERVICES 

O’  CY 

ATTN  5. 

BRECHT 

CO.,  INC. 

HMES 

BOEING  C 

QMPANY, 

THE 

COURT  STREET 

P.O.  30X 

3707 

SYRACUSE,  NY  13201 

j  i  A  .  .  L  s.  , 

W A  98! 

24 

0 1 C Y  ATTN  G.  MILLMAN 

01  CY 

ATTN  G. 

KEISTER 

0  1  CY 

ATTN  D. 

.  MURRAY 

01  CY 

ATTN  G. 

,  HALL 

01  CY 

ATTN  J, 

,  KENNEY 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS  ,  A.<  99’01 

:ALL  CLASS  ATTN:  SECURITY  OFFICER) 

0  1  C  Y  ATTN  T.N.  DAVIS  (JNCLASS  ONLY* 
D’CY  ATTN  TECHNICAL  LI3RARY 
0 1  C  Y  ATTN  NEAL  SHOWN  '.UNCLASS  ONLY) 

GTE  SYL/ANIA.  INC. 

ELECTRONICS  SYSTEMS  ORF-EASTERN  DIV 
7 T  A  STREET 
NEEDHAM,  MA  02194 

01 CY  ATTN  DICK  STSINHCF 


KAMAN  SCIENCES  CORP 
P  .0  .  3CX  ?463 

COLORADO  SFRINOS,  CC  3C932 
0  1  C  Y  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANC 
STUDIES 

312  STATE  STREET  P.C  DRAWER 
SANTA  3  A  R  3  A  R  A ,  CA  93102 
01 CY  ATTN  DAS  I AC 
Oi CY  ATTN  WARREN  3.  KNAP? 

0 1 C Y  ATTN  WILLIAM  MCNAMARA 
0 1 C  Y  ATTN  3.  GAM3ILL 


HSS.  INC. 

2  ALFRED  CIRCLE 
3EDF0RD,  MA  01730 

0 1 C  Y  ATTN  DONALD  HANSEN 


LINKABIT  CC RP 
10453  ROSELLE 
SAN  DIEGO,  CA  9212’ 

0 1 C  Y  ATTN  IRWIN  JACC3S 


ILLINOIS,  UNIVERSITY  OF 
107  COSLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN ,  IL  6 1 320 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 
0 1 C Y  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1901  NO.  3EAU REGARD  STREET 
ALEXANDRIA ,  VA  2231 1 
01 CY  ATTN  J.M.  A E IN 
01 CY  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  WOLFARD 
01 CY  ATTN  JOEL  3ENGST0N 


LOCKHEED  MISSILES  4  SPACE  CO. 
P.C.  BOX  504 
SUNNYVALE,  CA  94038 
0  1  C Y  ATTN  DEPT  50-12 
01 CY  ATTN  D.R .  CHURCH!  L(. 

LOCKHEED  MISSILES  4  SPACE  CO. 
3251  HANOVER  STREET 
PALC  ALTO,  CA  94304 

01  CY  ATTN  MARTIN  WALT  IE? 
01 CY  ATTN  W.L.  IMKIF  DEPT 
01 DY  ATTN  RICHARO  G .  JOHNS 
DEPT  52-12 

01 DY  ATTN  U.3.  CLAD  IS  IE 


INTL  TEL  4  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLS:.  NJ  0~’*0 

0  '  I Y  ATTN  TECHNICAL  LIBRARY 

J A  YCO  R 

11011  TORREYANA  ROAD 

P.O.  30 X  35154 

SAN  DI EGO  ,  CA  92133 

01 CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20310 

0 1 C Y  ATTN  DOCUMENT  LIBRARIAN 
01 CY  ATTN  THOMAS  POTSMRA 
01 CY  ATTN  JOHN  DASSOULAS 


MARTIN  MARIETTA  OCR? 
ORLANDO  DIV  IS!  IN 
P.O.  BOX  533' 
ORLANDO,  rL  3 ‘ : ? 


MCDONNSL  DOUGLAS  CORPORATION 
5301  30LSA  AVENUE 
HUNTINGTON  BEACH,  DA  925  -:' 
COY  ATTN  N.  HARRIS 
O’  CY  ATTN  J.  MOULS 
0  '  C Y  ATTN  GEORGE  MF.CZ 
O’CY  ATTN  W.  OLSON 
0 1 C  Y  ATTN  R.W.  HALPRIN 
0 ' C  Y  ATTN  TECHNICAL 

LI3RARY  SERVICES 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 


i  N  T  A 

BARBARA 

.,  CA  9310  1 

C  1  CY 

ATTN 

P.  FISCHER 

01  CY 

ATTN 

W.F.  C REV I  SR 

01  CY 

ATTN 

STEVEN  L.  GUT3CHE 

0  1  CY 

ATTN 

R.  30GUSCH 

01  CY 

ATTN 

R.  HENDRICK 

01  CY 

ATTN 

RALPH  KILB 

0  1  CY 

ATTN 

DAVE  SOWLE 

01  CY 

ATTN 

F.  FAJEN 

01  CY 

ATTN 

M.  SCHEI3E 

0  1  C  Y 

ATTN 

CONRAD  L.  LONGMIRE 

01  CY 

ATTN 

B.  WHITE 

01  CY 

ATTN 

R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE ,  NM  87106 
0 1  C  Y  R.  STELLINGWERF 
01  CY  M.  ALME 
0  1  C  Y  L.  WRIGHT 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 

1820  DOLLY  MADISON  BLVD 

MCLEAN,  V A  22101 

0  1C  Y  ATTN  W.  HALL 
O'CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 

1  2340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01 CY  ATTN  E.C.  FIELD,  JR. 


PHOTOMETRICS,  INC. 

•*  ARROW  DRIVE 
WOBURN,  MA  01801 

0 1  C Y  ATTN  IRVING  L.  K0F3XY 

PHYSICAL  DYNAMICS,  INC. 

?. 0.  BOX  3027 
BELLEVUE,  WA  98009 

01 CY  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94510 
ATTN  A.  THOMSON 

RID  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

0 1 C Y  ATTN  FORREST  GILMORE 
01 CY  ATTN  WILLIAM  B.  WRIGHT 
0 1 C Y  ATTN  WILLIAM  J.  K ARC AS 
01 CY  ATTN  H.  ORY 
0 1 C Y  ATTN  C.  MACDONALD 

RAND  CORPORATION,  THE 
1  5450  COHASSET  STREET 
VAN  NUYS,  CA  91406 

0 1 C Y  ATTN  CULLEN  CRAIN 
0 1 C Y  ATTN  ED  SEDROZI AN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01  CY  ATTN  3ARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  U2nd  STREET 
NEW  YORK,  NY  10036 

0 1 C Y  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS ,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

0 1 C Y  ATTN  LEWIS  M.  LINSON 
0 1 CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  S.  FRIEMAN 
0 1 C Y  ATTN  E . A .  STRAXER 
0^ CY  ATTN  CURT: 3  ,.  SMITH 

SCIENCE  APPLICATIONS  ,  INC 
1710  GOODRIDCE  DR. 

MCLEAN,  VA  22102 
0 1 C Y  J.  COCXAYNS 
01 CY  E.  HYMAN 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
3' 3  ELECTRICAL  ENGINEERING  EAST 
J N I  V  E  R  3  I T  Y  PARK,  ?  A  16802 
(NO  C-A33  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAB 


SRI  INTERNATIONAL 
323  RAVENSWOOD  AVENUE 


NLO  ? 

ARK. 

CA  9“025 

O’CY 

ATTN 

J.  CASPER 

0  1  0  Y 

ATTN 

DONALD  NSILSON 

2  i  -  v 

ATTN 

ALAN  3URNS 

O’CY 

ATTN 

G.  SMITH 

0  1  C  Y 

ATTN 

R.  TS’JNODA 

O’  CY 

ATTN 

DAVID  A.  JOHNSON 

0’  OY 

ATTN 

WALTER  G.  CHESNUT 

o’  :y 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J.  VICKREY 

01  CY 

ATTN 

RAY  L.  LEADABRAND 

01  CY 

ATTN 

G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

0  1  CY 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

01  CY 

ATTN 

D.  MCDANIEL 

TECH NO  LOOT  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
3EDF0RD,  MA  01T30 

O'CY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  i  SPACE  STS  GROUP 


ONE  SPACE  PARK 


REDONDO 

3EACH , 
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